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Generation, interaction and detection of dark solitons in 
Bose-Einstein condensates is considered. In particular, we 
focus on the dynamics resulting from phase imprinting and 
density engineering. The generation of soliton pairs as well as 
their interaction is also considered. Finally, motivated by the 
recent experimental results of Cornish et al. (Phys. Rev Lett. 
85, 1795, 2000), we analyze the stability of dark solitons un- 
der changes of the scattering length and thereby demonstrate 
a new way to detect them. Our theoretical and numerical 
results compare well with the existing experimental ones and 
provide guidance for future experiments. 



I. INTRODUCTION 

Bose-Einstein condensates (BEC's) offer a unique pos- 
sibility of studying nonlinear effects using matter waves. 
This has been spectacularly shown in the recent BEC ex- 
periments which demonstrate, among other things, the 
possibility of four wave mixing the creation of topo- 
logical structures such as vortices ||,^, the creation of 
solitons , as well as other demonstrations of the su- 
perfluid character 

Solitons are one dimensional waves that propagate 
without spreading in a nonlinear medium. Their shape 
remains unaltered after interacting with other solitons. 
The nonlinear Schrodinger equation, which accurately 
describes dilute BEC's at zero temperature, supports 
soliton solutions for attractive as well as for repulsive 
2-body interactions [||J^. These solutions correspond 
to macroscopically excited states of the mean field of 
the condensate. For a single component condensate, the 
3D time-dependent nonlinear Schrodinger equation, also 
known as the Gross-Pitaevskii equation (GPE), reads 
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Here g = ATifP'a/m, where a corresponds to the s-wave 
scattering length for binary collisions between atoms, m 
refers to the mass of the atoms and V to the trap po- 
tential. For repulsive interactions (a > 0), solitons are 
characterized by a local density minimum together with 



a sharp phase gradient of the wave function at the posi- 
tion of the minimum. In this case, the nonlinear effective 
mean field potential term g\^{r,t)\'^ balances the disper- 
sion of the wavefunction caused by the kinetic energy. 
Because there is a notch in the density these solutions 
are termed dark solitons. Only when the density of the 
condensate exactly vanishes at the density minimum, is 
the resulting soliton stationary. Such a standing soliton 
has exactly a phase jump of tt between the two parts of 
the condensate connected by it and its velocity is zero. In 
general, the local minimum density can range from max- 
imal to zero depth with the associated soliton velocity, q, 
ranging from zero to the speed of sound, Cs, in the con- 
densate, i.e. Q < q < Cs — \/ riQg/m, where rto refers to 
the mean density in the condensate. In a homogeneous 
ID condensate, a soliton solution can be analytically ob- 
tained in an elegant way by using the inverse scattering 
method |^. The wavefunction corresponding to a dark 
soliton located at q propagating along the z-axis with 
speed V is described by 
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In the case of a dark soliton propagating along the z-axis 
in an elongated 3D condensate, z refers to the position 
of the notch (x — y) plane, and v the constant velocity 
of such a plane with respect to a stationary background. 
The soliton size is of the order of twice the healing length 
lo = h/y/mgno. 

In contrast, for condensates with attractive interac- 
tions (a < 0) soliton solutions are characterized by a 
maximum in the density profile without any phase jump 
across it. These solitons are termed bright solitons, and 
the solution corresponding to a bright soliton reads: 
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state of a condensate with attractive inter- 
sufficiently large nonlinearity, i.e., g > gmim 
bright sohton Q. 

bright solitons described by the nonlinear 
equation have been extensively studied in 
of nonlinear optics (see and references 
condensates only dark solitons have very re- 
observed as nonlinear matter waves l0-p|. A 
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crucial difference between optical and matter wave soli- 
tons appears at first glance. While optical solitons are 
created in optical guides, i.e. in a cylindrical medium 
which is a priori unbounded, a condensate is always con- 
fined in a trap. To what extent the boundary conditions 
affect the properties and stability of the solitons has been 
the subject of recent theoretical studies (see e.g. ]Tl||). 

The literature concerning bright solitons in matter 
waves is not very extensive [p[p^-p^]. Condensates with 
attractive interactions in 2 and 3D are unstable objects 
that collapse very rapidly when the particle number be- 
comes too large [^,0. Therefore a bright soliton in a 
condensate with attractive interactions is also an unsta- 
ble object. However, since quasi- ID condensates with 
attractive interactions are stable, bright solitons could 
be generated in them JTst . Alternatively, an interesting 
approach based on vector solitons in two-component con- 
densates (with repulsive interactions) can also be used to 
to study bright solitons in 2 and 3D condensates. By cre- 
ating a dark soliton in one of the species one can induce 
a bright soliton-like structure in the other one |]6|, p^ , p^ 
which fills the minimum of the first species. Since a 
condensate with positive scattering length is stable, the 
bright soliton thus created is no longer limited by the 
collapse of the condensate. Finally, it has been recently 
shown that with the well-controlled use of Feshbach res- 
onances in *^Rb, it is possible to make a condensate with 
attractive interactions in a very controllable way . 
This opens a new way to study bright solitons in matter 
waves. 

The study of matter wave solitons, both experimen- 
tally and theoretically, principally involves three different 
aspects: their generation, coherent evolution including 
coherent effects during detection, and incoherent evolu- 
tion and dissipation. So far the observed solitons in one 
component EEC's have been generated by the method of 
phase imprinting. This method, originally proposed to 
generate vortices, has become a veryefRcient tool to en- 
gineer the phase in condensates [^2|j2^. Optimization of 
the phase imprinting method has been recently discussed 
by Carr et al. , where initially not only the phase but 
also the density is properly engineered. A proper combi- 
nation of both effects in a quasi- homogeneous condensate 
produces a stable, standing soliton whose properties can 
be used to test fundamental aspects of many-body theory 
such as quantum and thermal fluctuations. 

The coherent evolution of solitons refers to the evolu- 
tion before dissipation takes place. While a condensate 
has a lifetime on the order of seconds, the recent exper- 
imental results have shown that for solitons the lifetime 
is on the order of 15 milliseconds (in one component con- 
densates). For shorter times the Gross-Pitaevskii equa- 
tion provides an accurate description of the soliton dy- 
namics . In homogeneous systems and in the absence 
of dissipation and/or thermal phonons, solitons maintain 
a constant velocity. In elongated harmonic traps a dark 
soliton (and in 3D the nodal plane of the soliton) oscil- 



lates in the trap with a frequency Q. = ujgl \f2 where Wa 
corresponds to the axial trap frequency ||ll|. Coherent 
evolution also concerns the generation and interaction 
of pairs of solitons. This coherent soliton evolution has 
been recently studied in the context of two-component 
condensates, where new and rich dynamics appears [p^ . 

Dissipative effects include both dynamical and thermal 
instabilities. Dynamical instabilities are due to the fact 
that solitons are indeed one dimensional objects. When 
embedded in higher dimensions (i.e. 2 or 3D conden- 
sates) their stability strongly depends on the geometry 
of the condensate. For sufficiently elongated traps with 
a high transverse confinement it is not possible to excite 
the transverse modes of the trap and the soliton is then 
dynamically stable. On the contrary, for a looser trans- 
verse confinement the transverse modes can be excited, 
making the soliton plane bend and undergo a snake in- 
stability. As a result, the soliton decays into phonons 
and more stable structures such as vortices and vortex 
rings [|^,p5|,p6t , or even more exotic objects such as svor- 
tices ]271|. Depending on the trap geometry, dynamical 
instabilities may occur on a shorter time scale than ther- 
mal instabilities. One should point out here that the 
dissipatory behavior of solitons corresponding to dynam- 
ical instabilities can also be described, using the Gross- 
Pitaevskii equation, as coherent but unstable evolution 
11,11 ■ 

Thermal instabilities appear due to the fact that soli- 
tons are collective excited states of the mean field of the 
condensate, and, therefore, decay into the ground state 
within a finite time. The dissipation consists in scatter- 
ing of phonons on the soliton's notch plane. Studies of 
dissipation, which are related to the interaction of a soli- 
ton with a thermal cloud, demand one to go beyond the 
GPE and use the Bogoliubov-deGennes equations which 
describe collective modes such as phonons. Since a soli- 
ton can be regarded as a particle with negative mass, such 
scattering accelerates the soliton until it reaches the effec- 
tive sound velocity and vanishes. This scenario has been 
described by Fedichev et al. p^ ], and has been recently 
studied for 3D solitons in elongated traps by Muryshev 
et al. [ pO[ . Both theory and experiment indicate that the 
life time of the soliton due to thermodynamic instability 
is of the order of 15 ms. This contrasts with the lifetime 
of vortices which is of the order of a few seconds . 

This paper focuses on the generation and coherent evo- 
lution of dark solitons in one-component condensates in 
elongated traps. It is organized as follows: Section II 
addresses the issue of the generation of solitons. In Sec- 
tion III the problem of soliton detection is investigated 
and the dynamics concerning the opening of the trap are 
studied. Section IV is devoted to the interaction between 
solitons. We discuss therein the experimental conditions 
under which the effects due to interactions can be ob- 
served. In the latter, we study the stability of a soliton 
when the scattering length of the condensate is changed 
in a controllable way. This study is stimulated by the 
recent experiments in *^Rb by Cornish et al. pO| , pl| . Fi- 
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nally we present our conclusions in Section V. 



II. GENERATION OF DARK SOLITONS 

In this section we focus on the generation of dark 
solitons in one component Bose-Einstein condensates by 
three distinct methods: (a) phase imprinting, (b) density 
engineering, i.e., non-adiabatic changes of the potential 
confining the condensate and (c) combination of the pre- 
vious two methods. 

As already mentioned in the introduction, even for 
thermodynamically unstable states such as solitons the 
GPE provides an excellent tool with which to describe the 
coherent dynamics on the relevant timescale. To avoid 
the effect of dynamical instabilities we restrict our analy- 
sis to the following configurations: (a) cigar-shaped con- 
densates with a high aspect ratio and (b) quasi-lD con- 
densates. These quasi-lD condensates can be obtained 
for both harmonic and box- like potentials. For the for- 
mer, the radial confinement frequency is required to be 
much larger than the mean field interaction between par- 
ticles Q. For the latter, the healing length is required to 
be of the same order as the transverse box length. These 
quasi- ID condensates can be experimentally realized by 
loading a condensate from an elongated magnetic trap 
into a dipole trap created by a blue detuned Laguerre- 
Gaussian laser beam [pil,B2| . 



A. Phase Imprinting 

The method of phase imprinting consists of passing a 
short off-resonant laser pulse through an appropriately 
designed absorption plate and impinging it on a conden- 
sate. In this way one can imprint the desired phase struc- 
ture on the condensate and hence create a dark soliton 
11,11 ■ 

The Gross-Pitaevskii equation reduces to an effective 
one-dimensional nonlinear Schrodinger equation (NLSE) 
when the radial frequency is larger than the mean par- 
ticle interaction and when the longitudinal dimension of 
the confining potential is much longer than its transverse 
ones 

in-<^{z,t) = l- — —+V{z)+ 
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Here V{z, t) describes the interaction with the external 
laser, i.e., denotes the dipole potential generated by the 
far detuned laser pulse which acts only in one part of 
the condensate, and V{z) refers to the time-independent 
trapping potential which remains constant during the 
whole process. Let us review here how the phase imprint- 
ing can lead to the creation of a soliton . For a laser 
pulse duration shorter than the correlation time of the 



condensate Tcor = ^/fJ-, where /x is the chemical poten- 
tial, the wave function acquires a local phase factor e~"^ 
without changing the condensate's density profile. To 
generate the appropriate phase distribution which leads 
to a soliton, it is sufficient to use a potential which acts 
only on half of the condensate, e.g. 
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where f{t) is the temporal envelope of the laser pulse nor- 
malized to J f{t)dt=l. This potential imprints a phase 
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where Ig refers to the width of the potential edge, which in 
turn determines the steepness of the imprinted phase gra- 
dient at zq. Attainable experimental values correspond 
to a 10-90%-absorption width of the phase step. For this 
reason we use a factor of 0.45 in Eq.^. In accordance 
to this limit, the experimental values le correspond to 
> 2 /xm. Thus the phase of a dark soliton is composed 
of two areas of constant phase connected by a steep gra- 
dient. 

In an elongated cigar-shaped condensate, no matter 
how accurately the phase-imprinting method is imple- 
mented, it is not possible to make a standing soliton. 
The soliton thus generated will be a moving soliton whose 
speed and depth is directly related to and the ampli- 
tude of the imprinted phase Acj). Figure |^ shows the 
results of a numerical simulation for the time evolution 
of the density and phase of the condensate within the 
first millisecond after a phase with a phase gradient of 
A(j} = TT and /e = 2 /^m has been imprinted. 

The imprinted phase profile leads to a velocity field, 
Vz{z) = {h/m)d4)(z)/dz. During the evolution on a 
timescale of the correlation time, this velocity field leads 
to a reduction of density in the region dvz/dz > 
(z > Zq), whereas in the region dvz/dz < (z < zq) 
the density increases. After the minimum and the max- 
imum in the density have fully developed, they begin 
to back-react significantly, as may be seen in the evo- 
lution of the phase. The region of the phase gradient 
begins to change and leads to a change in the dynamics 
of the density distribution. The phase gradient splits up 
into two regions with phase gradients of similar shapes 
(A01 « A02 ~ A(/>/2). The density maximum is con- 
nected with one of the phase gradients and moves ap- 
proximately with the speed of sound, Cg = ^/Awnoah/m, 
towards negative z-values. As time increases the den- 
sity wave broadens due to dispersion and to the repulsive 
two-particle interactions. In contrast, for the minimum 
propagating towards positive z- values, the reduced inter- 
action energy results in a compensation of the dispersion. 
This leads to an increase in the steepness of the gradient 
together with a reduction of the width of the minimum. 
In this process, a second, less pronounced minimum to- 
gether with additional density perturbations are created. 
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As may be seen in Fig. the created pronounced min- 
imum, in connection with the tanh-phase distribution, 
propagates as a stable sohtary wave. 

The time scale needed for such a structure to develop is 
approximately given by Tds ~ Tcor ■ ('e/'o)j where Tcor and 
Iq correspond to the correlation time and healing length 
of the condensate, and the phase step connected with the 
soliton, A(f>2, accounts for approximately one half of the 
initially imprinted phase step A(f). 

The phase imprinting method depends very strongly 
on the width of the potential edge Zg as well as on the 
value of the imprinted phase difference A(j). For a width 
much larger than the healing length, l,, 3> ^O: the time 
needed for the dark soliton to arise is significantly en- 
larged, and only shallow solitons can be generated. For 
example, for a phase width of = 5 fim, but otherwise 
identical parameters to those used in Fig. (Q), the soli- 
ton structure develops only after an evolution time of 
tf^y « 15 ms. On the other hand, a phase imprinting with 
phases A$ > tt but with the same imprinting width l^, 
leads to a faster development of soliton structures, ac- 
companied always by the simultaneous creation of mul- 
tiple solitons. 



B. Creation of dark solitons by density engineering 

The possibility of creating solitons in a BEC by en- 
gineering only the phase suggests that it should also be 
possible to create solitons by purely engineering the den- 
sity distribution. This would be equivalent to the cre- 
ation of optical dark solitons by intensity modulations of 
a light field propagating in a nonlinear medium p3[ |. 

Experimentally it is simple to engineer the density of 
a condensate. For instance, one can modify the mag- 
netic trapping potential in which the BEC forms with an 
additional optical dipole potential of a far blue detuned 
laser beam, which is focused to form a thin "light-sheet" 
perpendicular to the long axis of the cigar-shaped con- 
densate. The spot size of the focus can easily be chosen 
to be much smaller than the axial size of the conden- 
sate. Therefore, in the Thomas-Fermi limit, the density 
distribution of the BEC will be described by an inverted 
parabola, except for the region where the laser focus is 
applied. This creates a local minimum, rimm, with a 
relative density, p = nmin/no-, which is controllable by 
the laser power. Here no is the density at the laser fo- 
cal position for negligible laser power, P = 0. For a 
Gaussian laser beam the shape of the density distribu- 
tion in the vicinity of the local minimum will be approx- 
imately an inverted Gaussian. The phase distribution 
remains constant over the whole condensate. In order 
to generate soliton structures the dipole potential is non- 
adiabatically switched off while the magnetic trap poten- 
tial is kept on. The phase and density distribution of the 
condensate adjusts to this new potential by creating pairs 
of equal but counter-propagating solitons. Note that the 



total phase over the condensate is conserved in this pro- 
cess. P should be ^ 0.01 in order to maintain phase 
coherence between the two portions of the condensate on 
either side of the density notch. 

We have numerically simulated the creation of dark 
solitons in a BEC in the range of parameters acces- 
sible to current experiments by pure density engineer- 
ing. For instance. Figure |(a) shows, 10 ms after the 
non-adiabatic removal of the optical potential, the den- 
sity and phase distribution of a ^^Na condensate with 
= 5 X 10^ atoms. In this case the magnetic trap has 
a radial trapping frequency ujj_ = 320 Hz and an aspect 
ratio of A = 25. On the other hand, the optical detuned 
laser has a Gaussian half width at 1 /e^ oiW — 2 ^m, and 
its intensity is assumed to be such that it corresponds to 
j3 = 0.9. The density profile of the condensate shows 
three pairs of counter-propagating dark solitons with dif- 
ferent velocities while the phase distribution depicts the 
corresponding steep phase gradients in the vicinity of the 
solitons. 

For a wider laser focus with respect to the heal- 
ing length, the number of dark soliton pairs increases; 
Fig. ||(b) shows the situation 10 ms after the switching 
off of a laser beam with (3 = 0.6 and W = 12 /im. One 
sees from the figure that after this time most of the initial 
density deformation has already been transformed into 
stable density minima. For a narrower laser focus with 
respect to the healing length it is possible to produce a 
single pair of solitons |24|. 



C. Phase and density engineering 

Finally, in this last subsection we briefly summarize the 
results which show that by a proper combination of phase 
and density engineering a single standing dark soliton can 
be created in a quasi-one dimensional BEC |Q. In this 
new scenario, a box like confining trap potential is used 
for the condensate. For a harmonic potential the Thomas 
Fermi radius scales as N^^^ and therefore the healing 
length scales as N~^^^, while for a box- like confinement 
the healing length scales as iV~^/^. Since the size of the 
soliton is of the order of twice the healing length it should 
be possible, for a box-like confinement with the appro- 
priate choice of experimental parameters, to dynamically 
observe the generated soliton in situ, i.e., without need- 
ing to first expand the condensate. The method of den- 
sity and phase engineering is simply a combination of the 
two above explained methods. First a density minimum 
is created by adiabatically ramping the intensity of a fo- 
cused laser into an initially uniform BEC (see Fig ^) . The 
focused laser field is abruptly switched off and a second 
far detuned laser pulse of uniform density is shined on 
one half of the condensate (phase imprinting). Thus the 
density minimum acquires the appropriate phase distri- 
bution. In this way one can create a single standing dark 
soliton. Variations on this technique allow one to create 
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in a well-controlled manner asymmetric soliton pairs and 
various combinations of larger numbers of solitons. 



III. OPENING OF THE TRAP 

In current experiments using phase imprinting the size 
of the created soliton ~ 2/o is, nevertheless, smaller than 
the diffraction limit of the wavelength of the imaging ra- 
diation, so the soliton cannot be observed in situ jj,^. To 
overcome this problem the trap is suddenly switched off, 
so that the condensate, and therefore the soliton, expands 
freely for a few milliseconds (tTOF) and thus becomes 
detectable via absorption imaging [Q. In this section we 
analyze the dynamics associated with the opening of the 
trap and the ballistic expansion. We consider a cigar- 
shaped geometry with a large aspect ratio. The ballis- 
tic expansion then occurs principally in the transverse 
direction. The dynamics related to the opening of the 
trap is complicated, since the abrupt switching off of the 
trap potential modifies not only the density distribution 
but also the phase structure present in the condensate. 
Our 3D numerical simulations show that as the conden- 
sate expands the soliton velocity diminishes very rapidly 
while its depth increases. Simultaneously a new mini- 
mum in the density distribution appears in the vicinity 
of the density maximum connected with one of the phase 
gradients. This new density minimum observed in the 
experiments Q travels opposite to the soliton direction 
with a velocity smaller than the sound velocity. Since a 
soliton can be interpreted as a particle with a negative 
mass, by opening the trap the soliton acquires kinetic en- 
ergy and, therefore, its velocity decreases until eventually 
it becomes a standing soliton. 

A condensate in a quasi-lD trap (infinitely long cylin- 
der along the axial direction) admits a scaling function 
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The corresponding stationary soliton solution reads 



v/r^tanh(zv/r^))e-*^''/* (8) 



where y — p/Rtf and z — z/Iq. We assume here for 
the radial coordinate that p, ~ Ptf- Switching off the 
trap abruptly corresponds to making ujp suddenly zero. 
A scaling solution takes the form 
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By substituting Eq. (||) into Eq. (|^) and splitting the real 
and imaginary part, a solution is found for an appropriate 
choice of the phase <?!), giving rise to the scaling 6z(t) — 

y^r~Hjj^<^. This approach, which is valid for uj^^ < 
tTOF ^ m/^'-^oj predicts a soliton velocity 



v{t) = v{0) 



ln(wpT + ^Jujjr^ + 1) 



(10) 



where t = t — topen so that w(0) corresponds to the veloc- 
ity of the soliton at the time the trap is suddenly switched 
off (topen)- This scaling law agrees very well with the nu- 
merical results obtained by solving the time dependent 
GPE, as is shown in Fig.(||), as well as with the experi- 
mental data of Ref. . 



IV. INTERACTING SOLITONS 

Solitons propagate in a nonlinear medium without 
changing their shape even when they interact with each 
other. The interaction between optical solitons has 
been intensively studied in the propagation of light in 
monomode optical fibers (see Ref. ||l^ and references 
therein) . Zakharov and Shabat demonstrated that homo- 
geneous Bose-Einstein condensates support multi-soliton 
solutions and that solitons interact with each other like 
classical particles with a short range repulsive interaction 
1^ . The signature of this repulsive interaction manifests 
itself as a negative shift in the position of each soliton 
after interaction as compared to the position of a single 
free moving soliton. For an untrapped homogeneous ID 
condensate this shift can be analytically calculated using 
the inverse scattering method. In this case the Gross- 
Pitaevskii equation reduces to 



1^2 

2'dz 



(11) 



where z is in units of the healing length and time in units 
of h/2gno- 

For a repulsive condensate, with the boundary con- 
ditions |5'(z,f)p ^constant, a soliton solution moving 
with constant velocity through the condensate can be re- 
expressed as [| 36 



^{z,t) 



(A 



■ exp(2i/(z — zq ~ Xt)) 



1 + exp(2j/(z — zq — Xt)) 



(12) 



where + = 1 and zo is the initial position of the 
soliton at i = 0. The parameter A characterizes the am- 
plitude and the velocity of the soliton in units of Cg. In 
these units — 1 < A < 1, where A = ±1 corresponds to 
a completely filled soliton (zero depth) moving with the 
speed of sound, whereas A = corresponds to a standing 
soliton. 
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To calculate the spatial shift due to a soliton-soliton 
interaction one simply compares each final soliton posi- 
tion to what it would have been had it not undergone a 
collision. For two solitons with velocities Ai and A2, the 
resulting shifts are given by P,p6| 



Szi 



SZ2 



1 
1 

1^2 



log 



■log 



(Ai - Az)^ + jvi + 1^2^ 
(Ai-A2)2 + (1.1-1.2)2 J 
(Ai - A2)2 + {lyi + 1^2) 
[Xi - X2r + [i^i - 1^2^ \ 



If the solitons have equal velocities, i.e. Ai = — A2, 
shift is the same for both of them: 



Sz = 



log[|Ai 



(13) 

(14) 
the 

(15) 



Generally speaking, to see how the inhomogeneity due 
to the trapping potential affects the interaction dynamics 
between solitons one has to turn to numerical solutions of 
the Gross-Pitaevskii equation. To this end, we solve the 
GPE when two opposite phase gradients are imprinted 
in a cigar-shaped condensate 



(z) — [tanh((z — zi)/l^i) — tanh((z - 



Z2)/le2)], (16) 



where Zi and l^i denote the positions of the phase gra- 
dients and the width of the potential edge respectively. 
Thus, at i = the wavcfunction of the condensate reads 



«'(z,< = 0) = 



V(^) 



(17) 



where V'(^) is the ground state solution of the Gross- 
Pitaevskii equation. After a time of the order of the 
correlation time, such a phase distribution generates two 
counter propagating solitons (two notch planes located 
at Zi and Z2) moving with velocities (Ai = — A2), to- 
gether with two counter propagating density waves (den- 
sity maxima) that move with the speed of sound. The 
soliton positions are monitored by following the density 
notch in the condensate. Figs. || and ^ display the results 
of a full 3D calculation corresponding to the experimen- 
tal conditions discussed in where the soliton planes 
are initially separated by 20 and 35 microns respectively. 
Due to the change in the background density of the con- 
densate, the solitons first accelerate, and then cross each 
other. In Fig. ^ for the case of an initial separation of 
20 microns the two solitons overlap strongly during the 
collision, which lasts for approximately 1 ms. In Fig. |^ 
for an initial separation of 35 microns, the solitons collide 
with such a high velocity, that the short interaction time 
makes the position shifts negligible. 

In spite of the 3D character of the elongated trap one 
can still use the analytical results ^ from the homoge- 
neous ID case by applying a local density approxima- 
tion where the condensate density is considered constant 
in the region of the soliton plane. This approximation 
is valid away from the the edges of the trap. Insert- 
ing the parameters used in the Hannover soliton experi- 
ment H for s^Rb with N=10^ a=5.7 nm, and an initial 



Vsoiiton = S.3^m/ms, and assuming an effective cross 
section area of S' = 25/1771^, the predicted shift in the 
position of the solitons is Sz — —O.lfim. This calcu- 
lated value is in agreement with the numerical simula- 
tions. Such a small positional shift cannot be detected 
experimentally. In order to obtain a shift of the order of 
10/im the velocity difference should be — W2I ~ 10~^Cs. 

Finally, in FigJ| we display the results (ID) for the 
collision between two dark solitons created with phase 
gradient of 1^ = 10. The initial velocities are very 
small (w ^ O.OScs). Thus the interaction time becomes 
very large. The solitons are clearly seen to bounce off 
each other like classical particles undergoing an elastic 
collision. The shift in the positions is, however, still 
very small. Here we have chosen the density to be 
4 X lO^^cm"'^ which again corresponds to a condensate 
with the effective cross sectional area S = 25fj,m^. 

In order to unambiguously detect experimentally a sig- 
nature of the interaction between solitons in condensates, 
one should create two solitons with very similar veloci- 
ties propagating in the same direction. The solitons will 
thus interact for a long time. From Eqs. ( p^ and ( p^ ) 
it is clear that this scenario can amount to an arbitrar- 
ily large shift. This requires, however, a very long con- 
densate. Such condensates are becoming experimentally 
available |3^,^^ in reduced geometries. On the other 
hand, this drawback could be removed by creating the 
condensate in a ring geometry. Another possibility is to 
create solitons in a hard wall geometry, as for example in 
a hollow blue-detuned laser beam with laser light sheet 
endcaps, so that the solitons will reflect from the ends 
without changing their form, since the wavcfunction at 
the endcaps acts as a pinned soliton |38|. 



V. STABILITY OF DARK SOLITONS FOR 
NEGATIVE SCATTERING LENGTHS 

Recent experiments concerning the use of Feshbach res- 
onances to change both the magnitude and the sign 
of the scattering length of condensates with alkali atoms 
offers a new range of phenomena to study. In particular, 
a carefully controlled study of the outstanding problem of 
collapse of the condensate for negative scattering length 
becomes possible As already mentioned in the pre- 
vious sections, the form of a soliton depends on the sign 
of the scattering length. For positive scattering lengths 
the stable soliton solution is a density notch, i.e, a den- 
sity minimum; for sufficiently strong negative scattering 
lengths the stable solution is a bright soliton which is a 
density peak. 

The possibility of changing the scattering length from 
positive to negative values opens the question of the sta- 
bility of dark solitons in attractive condensates. Let us 
first reconsider the stability of dark solitons when adia- 
batically changing the scattering length. Since a soliton 
is a particular solution of the Gross-Pitaevskii equation 
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for a well defined scattering length, it is reasonable to 
assume an adiabatic change of the scattering length will 
change the velocity and depth of the soliton gradually. 

On the contrary, an abrupt change in the condensate's 
scattering length will destroy the soliton. For negative 
scattering lengths the physical and mathematical situa- 
tion changes dramatically. If we only consider low den- 
sities and neglect three-body recombination - which be- 
comes important at high densities and produces addi- 
tional kinetic energy - in the absence of any solitons, 
the instability of the condensate is seen as a collapse of 
the condensate's wavefunction. This is shown in Fig. 
^, where one can see a shrinking cloud as time pro- 
ceeds and where no solitons are present. In these sim- 
ulations we have used a small ®^Rb condensate with 
5000 atoms in the same cigar-shaped trap previously dis- 
cussed. We change the scattering length from its initial 
value a = 5.7nm to a —> —0.1a. In the presence of a dark 
soliton, (Fig. the scenario changes dramatically. The 
soliton splits the cloud into two separate parts which in- 
dependently continue to collapse. A direct consequence 
of non-adiabatically changing the sign of the scattering 
length in the presence of a soliton is the creation of a 
large number of density waves. This effect speeds up the 
collapse of the wave function because of the local increase 
in the density. A more careful study which takes into ac- 
count three-body recombination processes is needed to 
investigate the dynamics for longer times. 



VI. CONCLUSIONS 

We have discussed the generation, evolution, and inter- 
action of dark solitons in matter waves. We have first re- 
viewed different approaches to generate standing or mov- 
ing dark solitons in one component condensates. The in- 
teraction dynamics between dark solitons have also been 
addressed and we have discussed under which circum- 
stances the interaction can be observed experimentally. 
We conclude that in present experiments using cigar- 
shaped condensates with a large aspect ratio, a conclusive 
signature of the soliton interaction cannot be observed. 
However, by using other geometries, such as quasi-lD or 
toroidal condensates, the interaction could be unambigu- 
ously detected experimentally. A stationary wave in the 
form of a density notch or peak, even if it moves with 
less than the speed of sound, cannot truly be called a 
soliton until it is demonstrated that it interacts as one, 
since that is the defining characteristic which gives soli- 
tons their particle-like nature. 

Finally, we have discussed the stability of dark solitons 
to sudden changes of the sign and value of the scattering 
length. We find that the presence of a dark soliton can be 
unambiguously detected by the radical change in the dy- 
namics of a collapsing cloud when a Feshbach resonance 
is used to tune the scattering length negative: the cloud 
splits in two. 
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FIG. 1. Numerical simulations (ID) showing the time 
evolution of the density profile of a ^'^Rb condensate with 
TV — 10^. The static trap has a frequency ujz ~ 2-k x 14 
Hz. The density profile n{z) and the phase are depicted re- 
spectively in (a) and (b) within the first millisecond after a 
phase-imprint with A(j) = vr and le = 2 /im has been imple- 
mented. 



FIG. 2. Time evolution of the density profile n{z) and 
the phase (j}{z) for the first 5 milliseconds after a 7r-phase has 
been imprinted (otherwise the same parameters as in Fig hi). 



FIG. 3. Evolution of pairs of counter propagating dark 
solitons created by a non-adiabatic change of the BEC 
trapping potential. The parameters here correspond to a 
Na-condensate with N = 5 x 10^ atoms in a magnetic trap 
of the clover-leaf type with a radial trapping frequency of 
uj± = 320 Hz and an aspect ratio of A = 25: (a) W = 2 /im, 
/3 = 0.9 and (b) /3 = 0.6 and W = 12 /im. In the same figure, 
density profile and phase gradients are depicted. Note that 
the phase in (b) has been plotted modulo 27r. 



FIG. 4. (a) A combination of a box-like potential and a 
tightly focused, blue-detuned laser beam is used to engineer 
the density, (b) The resulting wavefunction is phase engi- 
neered with a second, far-detuned laser beam, resulting in an 
initial state that resembles very much a standing dark soliton. 



FIG. 5. Velocity of the soliton outside the trap versus ob- 
servation time. The soliton is imprinted at t = and the trap 
is suddenly removed after t — 4 ms. The curve shows the an- 
alytical scaling law, and the dots correspond to the values of 
the soliton velocity obtained from numerical 3D simulations. 



FIG. 6. Interaction between solitons in a cigar shaped con- 
densate. The solitons are initially created by the method 
of phase imprinting with an initial separation of 20 microns. 
Here N=10^ and a=5.7 nm. 



FIG. 7. The same parameters as in Fig. 6, with the solitons 
initially separated by 35 microns. At the time of the collision 
their velocities are now higher, hence the interaction time 
shorter. 



FIG. 8. The interaction between solitons in the ID case. 
The density waves created by the phase-imprinting are mov- 
ing with the speed of sound and are reflected at the condensate 
boundary (not shown in the figure) which eventually results 
in crossing waves at t — 12ms. 
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FIG. 9. The scattering length is changed at t = 2ms from 
a to —0.1a. The cloud, represented here as the integrated 
density as a function of z, is starting to collapse at the onset 
of the negative scattering length. 

FIG. 10. Same situation as in Fig. ^ with the presence of a 
dark soliton. The soliton splits the cloud into two parts which 
independently start to collapse. 
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